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MICHAEL BAAKE 

Abstract. Discrete point sets S such as lattices or quasiperiodic Delone sets may permit, 
beyond their symmetries, certain isometries R such that S n RS is a subset of <S of finite 
density. These are the so-called coincidence isometries. They are important in understanding 
and classifying grain boundaries and twins in crystals and quasicrystals. It is the purpose 
of this contribution to introduce the corresponding coincidence problem in a mathematical 
setting and to demonstrate how it can be solved algebraically in dimensions 2, 3 and 4. 
Various examples both from crystals and quasicrystals are treated explicitly, in particular 
(hyper-)cubic lattices and quasicrystals with non-crystallographic point groups of type H2, 
H3 and Hi. We derive parametrizations of all linear coincidence isometries, determine the 
corresponding coincidence index (the reciprocal of the density of coinciding points, also 
called S- factor), and finally encapsulate their statistics in suitable Dirichlet series generating 
functions. 



1. Introduction 

The concept of a coincidence site lattice (CSL) arises in the crystallography of grain and 
twin boundaries |151 1431 138] . Different domains of a crystal across a boundary are related 
by having a sublattice (of full rank) in common. This is the CSL. It can be viewed as the 
intersection of a lattice with a rotated copy of itself, where the points in common form a 
sublattice of finite index (we shall not discuss any situation other than that). Until recently, 
CSLs have been investigated only for true lattices or for crystallographic packings, for example 
cubic or hexagonal crystals (SIEHIEZI- Although the subject itself is quite old, no systematic 
investigation in more mathematical terms has been carried out so far. 

With the advent of quasicrystals, many new cases arose, since quasicrystals also have grain 
boundaries and one would like to know the coincidence site quasilattices HH1 E3 02] and, 
more specifically, which of them can form twins (or multiple twins, where the angle between 
the grains is a rational multiple of ir). Added impetus is given by the experimental progress 
made in recent years [20], in particular on the study of twins in icosahedral structures. So, 
an extension of the CSL analysis to all discrete structures is desirable. 

It is the aim of this article to provide the mathematical basis for it, and to display various 
relevant examples in detail. Though we shall include many proofs, it is not possible to present 
a complete account here, and we have to refer to original sources several times. Additional 
material can be found in ^2 EH] IB] • Let us briefly outline how this article is organized. 
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In Section 2, we set the scene for the periodic case and derive various results on the 
coincidence structure of lattices. As we shall need a generalization to modules later on, most 
results are formulated to match that. 

Section 3 deals with cubic lattices, where the cases of dimensions 2, 3 and 4 are treated 
explicitly. The groups of coincidence isometries are derived together with the corresponding 
index formulas, and the statistics of the CSLs is encapsulated in Dirichlet series generating 
functions which are related to various Dedekind zeta functions. 

Sections 4 and 5 consider quasicrystals. Here, the problem must be split into two parts, one 
being the coincidence problem for the underlying limit translation module (which is universal 
and discussed in detail) and the other being a problem specific correction in the projection 
formalism (which is only briefly outlined). This is followed by some concluding remarks. 

Finally, in the Appendix, a closely related problem is presented. While CSLs are special 
sublattices which depend on metric properties, the number of all sublattices of a given index 
is an affine property and depends only on the rank. The solution is derived explicitly for 
arbitrary rank. 

2. Preliminaries and some general results 

This paragraph focuses on lattices, although many properties could directly be formulated 
for modules. We prefer this approach, as not every reader might directly want to go beyond the 
lattice situation. We will then generalize the concepts briefly when we pass on to quasicrystals. 

The first concept we need is that of a lattice in d-dimensional Euclidean space (or d-space 
for short), where we follow the standard approach, see JJj or [2Hj f° r details. 

Definition 2.1. A discrete subset r of E d is called a lattice (of rank and dimension d) if it is 
the Z-span of d vectors Oj, . . . , a d that are linearly independent over M. These vectors form 
a basis of the lattice. 

In particular, we can write r = 7La x © • • • © Za d , and r is isomorphic to the free Abelian 
group of rank d. Another rather common (and equivalent) characterization is to say that r 
is a co-compact discrete subgroup of M d . 

Beyond a lattice r, we shall also need its dual, r*, which is given by 

Definition 2.2. T* := {x | xy G Z for all y G T}. 

Here, x y denotes the standard Euclidean scalar product. The dual is of course also a 
lattice, since a\,...,a* d (defined through a* k ■ a £ = Ski) 1S a basis of r*, called the dual basis. 
It is convenient to attach a basis matrix B to a lattice r, where the fc-th column of B consist of 
the coordinates of a& in the standard Euclidean basis e 1; . . . , e d . It is clear from Definition 12. 21 
that, if B is a basis matrix for r, then B* := (B~ 1 ) t is one for r*. 

To proceed, we need the concept of a sublattice. Here, in view of generalizations to come 
in later chapters, we employ the group structure. 

Definition 2.3. Let r be a lattice in K d . A subset r' C r is called a sublattice of r when 
it is a subgroup of finite index. The latter is the number of residue classes (or cosets) of r' 
in r, denoted by [r : 
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By definition, a sublattice has finite index in its (super-)lattice, and thus shares rank and 
dimension with it. In particular, Z 2 is not a sublattice of Z 3 in this notation. If such a 
situation occurs, we shall speak of net planes or lattice planes (as Z 2 = Z 3 n {z = 0}). 

For lattices, the index also has a direct geometric meaning: it is nothing but the (inverse) 
quotient of the volumes of the fundamental domains of the two lattices, as can easily be seen 
with explicit bases for them. A well known result |17j in this context is the following. 

Lemma 2.1. Let r be a lattice in K d with basis matrix B. Then, T' is a sublattice of T if 
and only if there exists a non-singular integral matrix Z such that B' = BZ is a basis matrix 
for r' . The corresponding index is [T : T'] = |det(Z)|. □ 

This description of sublattices is indeed helpful as one can now show 

Lemma 2.2. Let i~" be a sublattice of T C E d of index m. Then, mF is a sublattice of T' of 
index m d ~ l . 

Proof: Note that T' C T means B' = BZ with Z integral and |det(Z)| = m, by Lemma ITT1 
But mB = mB'Z^ 1 = B'Z' is a basis matrix for mT, and Z' = mZ~ 1 is integral (by the 
standard formula for the inverse of a matrix). This means mT C i~" by Lemma 12.11 Finally, 

gives the statement about the indices. □ 

Note that, for d > 1, mT need not be the maximal sublattice of T' which is a homothetic 
copy of T, as can be seen from the example T = Z 2 , T' = 2Z 2 . Concerning the dual lattice, 
one can show another useful result. 

Lemma 2.3. If T<i is a sublattice of T\, one has [T£ : r*] = [Ji : Fz], and the corresponding 
factor groups are isomorphic: T^/T^ ~ T 1 /r 2 . 

Proof: By assumption, B<i = B\Z with Z integral and non-singular. Consequently, one has 
B\ = B 2 Z t and the first statement follows from det(Z*) = det(Z). The two Abelian factor 
groups are thus of equal order. Isomorphism follows from the observation that the subgroups 
of i~i which contain i~2 are, by duality and part one of the Lemma, in one-to-one relation with 
the subgroups of T 2 which contain T^, with preservation of the corresponding indices. □ 

Having prepared the ground, we can now introduce various concepts and results which will 
be helpful for the coincidence problem. 

Definition 2.4. Two lattices ri,i~2 are called commensurate, denoted by Ti ~ I~2, when 
r% n i~2 is a sublattice (of finite index) of both T± and i~2- 

It is clear that only lattices of the same rank can be commensurate to one another, and 
one also has 

Proposition 2.1. Commensurateness of lattices is an equivalence relation. 

Proof: Reflexivity and symmetry are clear by definition. Finally, T\ ~ i~2 and i~2 ~ ^3 
together imply that T\ fl J2 H I3 is of finite index in both T\ and Z3 (one way to see this is 
via suitable use of Lemma I2.2J) , which implies T\ ~ 7~3 and hence transitivity. □ 
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Proposition 2.2. Let A and A ^ e commensurate lattices. Then, A n A an <^ A + A are 
lattices, and one has the following diagram, 

A + A 
A r 2 

«1 \ / n 2 

An A 

where A <Z B is written as A — > 1? and i/ie indices satisfy m 1 = n 2 and m 2 = n l . 
Furthermore, the following equations hold: 

(A n A)* = A* + A* and (A + A)* = A* n A* • 

Proof: Recall that A + A is defined as 

A + A := {x x +x 2 \x l er i , x 2 e A} • 

It is the smallest group which contains the lattices A and A • As A ~ A , Lemma 12.21 
guarantees that A + A is again a lattice because fc(A + A) C A H A for some k £ N. 

The claim about the indices is a direct consequence of the second isomorphism theorem for 
groups giving 

A/(A n A) =i (A+AVA and A/(A n A) ^ (A + AVA- 

As A H A is of finite index both in A and A> all these factor groups are finite. 

The second claim follows directly from the definition of the dual lattice: x £ (A + A)* 
means x ■ y S Z for all y £ (A + A)i hence in particular for all y G Ai & = L 2, and therefore 
a; G (A* n A*)- Conversely, if x ■ y k G Z for y fc £ T k , k £ {1, 2}, then a; ■ + y 2 ) G Z and 
(A n A*) C (A + A)*- The other identity follows by duality. □ 

Now we come to the central concept of this article. 

Definition 2.5. Let r be a lattice in E . An orthogonal transformation R £ 0(d) := 0(d, M) 
is called a coincidence isometry of r when Rr ~ r. The integer £(R) := [f : (r n -RA^] i s 
called the coincidence index of R. If i? is not a coincidence isometry, X'(-R) := oo. We set 

OC(r) := {i? G 0(d) I < oo} 

soc(r) : = {r £ oc(r) | det(i?) = 1} 

If necessary, we shall use the lattice as a subscript for the coincidence index, U = Up, but 
usually we can safely manage without. In this article, we consider only linear isometries, called 
isometries for simplicity. In more general situations, extensions to affine transformations are 
necessary, compare [313 App. A], but we shall not use them here. One immediate consequence 
of Definition 12.51 is 

Theorem 2.1. If T is a lattice in E d , OC(A) and SOC(A) are subgroups of 0(d). 
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PROOF: Let R± and R 2 be in OC(T), i.e., r ~ R±r and r ~ i^-f- Clearly, this also implies 
r ~ R 2 1 T and ~ R^R^ T. From transitivity, we may conclude that r ~ R l R 2 ~ 1 r and 
hence R X R^ G OC(r). So, OC(r) is a subgroup of 0(d); the corresponding statement for 
SOC(r) is obvious. □ 

Without further information on T, one cannot determine the corresponding index £ = Up, 
except that we know: 

Lemma 2.4. Let R G OC(T). Then, R' 1 G OC(r) and A7(i?) = S(R^ 1 ). 

Proof: The first claim is clear from the group property of OC(-T). Since R is an isometry, 
its action does not change the volume of fundamental domains, and we obtain the equation 

e(r) = [r-.(rnRr)] = [Rr-.(mRr)} = [r ■. (R^r n r)] = ^(ir 1 ). 

This establishes the claim. □ 
More can be said about the OC-groups of related lattices. The following is immediate: 

Lemma 2.5. If T is a lattice, R an orthogonal transformation and A G M\{0}, one has the 
relations OC(\r) = OC(r) and OC(Rr) = ROC^R- 1 ~ OC(r). □ 

Since orthogonal transformations and homotheties commute, the OC-group does not change 
(up to conjugation) if one applies linear similarity transformations. This is one step to show 
that the OC-group essentially is an invariant of so-called Bravais classes (cf. [3U] for details on 
this concept) - if one disregards non-generic solutions for special representatives. The latter 
problem does not occur if one deals with irreducible symmetries, such as fourfold symmetry 
in the plane or full cubic symmetry in E 3 . In this article, no other situation shall be discussed. 

Lemma 2.6. Let Fi be a sublattice of F\ of (finite) index m = [Ji : i^]. Then, 0C(i~2) = 
OC(-Ti) and, for any isometry R out of this group, one has £\(R) \ mX^i?). 

Proof: If R is a coincidence isometry of i~2> we know that i~2 H RJ2 has finite index in i~2 
and hence also in T 1 . As F 2 n RF 2 C F x D RF X C A, we see OC(r 2 ) C OC(A). On the other 
hand, \T\ : F 2 \ = m implies mF\ C F 2 by Lemma so also OC(i~i) C OC(r 2 ) and the two 
groups must be equal. The second statement follows from T 2 n RT 2 C T% fl i?A C A and 
i~2 H ii/2 C i~2 C i"i by a comparison of indices. □ 

If we write T\ = lj^=o + ^2) with to = 0, one can say more about the relation between 
S\(R) and £ 2 (R) if R respects the cosets in suitable ways. In particular, if RT 2 is disjoint 
from the cosets ti + F 2 for all i > 0, one can also derive that U 2 (R)\Ui(R). 

The last results on the coincidence groups can obviously be extended to 

Corollary 2.1. Commensurate lattices possess the same OC-group. □ 

With the above results, one can now also relate a lattice with its dual. 

Theorem 2.2. Let F be a lattice in E d and T* its dual. Then, OC(r*) = OC(r) and the 
coincidence index of any orthogonal matrix is the same for both lattices. 
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Proof: Assume F ~ FF. Then, since F** = F and (Rr)* = Rr*, we have oo > m = 
27(12) = [F : (F n FF)] = [(r + AT) : F] = [F* : (F + Rr)*} = [F* : (F* n 12F*)] = 27*(F), 
by applying Proposition 12.21 (specialized to the case F2 = 12Fi, where m 1 = m 2 = n l = n 2 ) 
and Lemma 12.31 It also follows that 27(F) = 00 if and only if 27* (F) = 00. □ 

We have now all prerequisites to tackle the coincidence problem for lattices and crystals. 
While we proceed, we shall need, in particular in the part on quasicrystals, various results 
from algebraic number theory. It is not possible to present a self-contained description here, 
but references to the relevant literature shall be given. For general background, we refer to 
[HJ Ch. 1.4], to [SOI Chs. XIV-XVII and Ch. XX], and, for a number theoretic approach to 
quasicrystals in general, to [39] . 

3. Lattices and crystals: the cubic case 

Although there are several cubic lattices, let us first consider the primitive (hyper-)cubic 
lattice, Z d . With the standard Euclidean basis e l3 . . . , e it can be written as 

Z d = Ze 1 e---®Ze d . 

As we shall need the term 'primitive' later on in a different meaning, we replace it here by 
the term F-type from now on, as compared to F-type (for face centred cubic, fee, and its 
generalizations) and to F-type (for body centred, bee, and its generalizations). 

It is immediately clear that an orthogonal matrix R with rational entries only is a coinci- 
dence isometry of Z d : one can directly give a lattice which is a common sublattice to both 
Z d and RZ d , namely mL d , where m is the denominator of R defined through 

(3.1) den(F) := gcd{k \ k ■ R integral} , 

where gcd is the greatest common divisor. On the other hand, as soon as one entry of R is 
irrational, Rij say, no point in the direction of ej coincides with a lattice point after rotation. 
This means that we have one lattice direction without any coincidence, hence infinitely many 
residue classes of the set of coinciding points which therefore is not of finite density, and R 
cannot be a coincidence isometry. So we have 

Theorem 3.1. OC(Z d ) = 0(d,Q) and SOC(Z d ) = SO(d,Q). □ 

Having this result for the F-type (hyper-)cubic lattices, the obvious next question is what 
happens for the other lattices with (hyper-)cubic symmetry. The types of such lattices, up to 
similarity transformations, are as follows. 

Theorem 3.2. In dimensions d = 3 and d > 5, there are precisely three (hyper-) cubic 
lattices, namely F-type, P-type and B-type. For d = 1 and d = 2, there is only one such 
lattice (represented by the integers resp. by the square lattice), while in d = 4, there are two 
such lattices (P-type hypercubic and centred). □ 

A proof can be found in (HOI- It is possible, in any dimension, to realize the (hyper-) 
cubic lattices in such a way that they are commensurate - usually, one works with as 
representative of F-type, with the root (weight) lattice D d (D d ) as representatives of F-type 
(F-type) in dimensions d > 5, with F 4 as centred lattice in 4-space, and with the usual fee and 
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bcc lattices in 3-space (or A 3 resp. A%), compare ^Sj. With this convention, the OC-groups of 



might differ. At present, we do not know the general answer in higher dimensions, but more 
can be said about the cubic lattices in dimensions 2, 3 and 4. So, let us summarize some of 
those results. 

3.1. d = 2: the square lattice. Let us first describe the case of the square lattice in more 
detail, as this is the simplest non-trivial example. Here, we shall be less formal and refer 
for all the details and proofs to [101 without any further mentioning. First, we restrict the 
description to rotations, and come to reflections at the end of the section. 

The square lattice Z 2 , embedded in E 2 resp. M 2 , consists of all integer linear combinations 
of the two vectors e± and i.e., Z 2 = r Le 1 © Ze 2 . A rotated copy Kl? with R £ SO(2,M) 
results in a CSL of finite index, according to Theorem 13.11 if and only if both cos^) and 
sin(</?) are rational, where (p is the rotation angle. This gives the well-known relation between 
coincidence rotations and primitive Pythagorean triples [31]. The group of coincidence rota- 
tions is thus explicitly seen to be SOC(Z 2 ) = SO(2,Q). To investigate this group, we employ 
some elementary results from the algebraic theory of quadratic fields 30 , Chs. XIV and XV]. 
In particular, we notice that, with i = \/— 1, we can identify Z 2 with the ring of Gaussian 
integers (the algebraic integers of the quadratic field an extension of Q of degree 2): 



The ring Z[i] is a Euclidean domain and thus has unique factorization up to units, see |3Uj . 
The units (i and its powers) form a group isomorphic to C4, the rotation part of Aut(Z 2 ). 

In this setting, a rotation R((p) £ SOC(Z 2 ) corresponds to multiplication by a complex 
number e lip G Q(i). This number can be written as e tip = a/ (3 with a,/3 £ Z[i] coprime and 
of equal norm, i.e., |a| = \/3\ (so, e lip rotates the lattice point (3 into another lattice point, 
a, on the same circle around the origin). We now factorize numerator and denominator into 
Gaussian primes. They are the rational (or ordinary) primes p = 3 mod 4, the factor 1 + i of 
2 (a so-called ramified prime), and the pairs of complex conjugate factors of rational primes 
p = 1 mod 4 (where p = LOptJ p , and uj p /uj p is not a unit in Z[i]). Clearly, as a and f3 are 
coprime and both divide the same rational integer i = \a\ 2 = |/3| 2 , only the last type of primes 
can occur in the factorization, always one (or a power of it) in the numerator and its complex 
conjugate in the denominator. 

As a consequence of the unique factorization property up to units, every coincidence rota- 
tion can then be factorized as 



where n p S Z (only finitely many of them ^ 0), e is a unit in Z[i] (a power of i), p runs 
through the rational (or ordinary) primes congruent to 1 (mod 4), and the iv p , U p are the 
(complex conjugate) Gaussian prime factors of p. We thus have the (non-trivial!) result that 
SOC(Z 2 ) is an infinitely generated Abelian group that nevertheless permits the factorization 



the (hyper-)cubic lattices in K d 



(with d fixed) are identical, though the corresponding indices 



(3.2) 



1? = Z[i] = {m + ni\m,n G Z} . 



(3.3) 
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into a torsion group and a free Abelian group, namely 

(3.4) SOC(Z 2 ) = S0(2,Q) ~ C 4 x 

with generators i for C4 and lo p /uj p with p = 1 (4) for the infinite cyclic groups. By Z^ ^ we 
mean, as usual, the infinite Abelian group that consists of all finite integer linear combinations 
in the (countably many) generators. We prefer a multiplicative rather than additive notation 
here as the coincidence groups will not be Abelian in later examples. 

To determine the coincidence index £(R), we observe that it is 1 for the units (i.e., true 
symmetry rotations) and p for the generator LO p /uJ p (since p = lo p lo p = norm(w p ) counts 
the number of residue classes of the corresponding CSL in Z 2 ). More generally, due to 
multiplicativity, 2J(R) is the (number theoretic) norm of the numerator of ()3.3|) . i.e., 

(3.5) E(R) = 11 p M . 

p=l (4) 

This shows the power of the generator approach in this (Abelian) situation. The first few 
generators with U > 1 are 

4 + 3i 12 + 5i 15 + 8z 21 + 20* 35 + 12i 40 + 9i 
5 ' 13 ' 17 ' 29 ' 37 ' 41 ' 6tC ' 
These are normalized (by multiplication with a suitable unit) to have argument in (0,7r/4), 
and are shown with denominator £ (a prime = 1 (mod 4)). All other coincidence rotations 
are obtained by combinations, as indicated in Eq. (|3.3I) . 

It is convenient to summarize the possible coincidence indices and the number of CSLs with 
a given index by means of a generating function. To do so, let 4/(m) denote the number of 
coincidence rotations of index m, which means that f(m) counts the different CSLs of index m. 
As a consequence of unique factorization in Z[i], /(m) is a multiplicative arithmetic function 
(i.e., /(I) = 1 and /(mira2) = f (mi) j (7772) for coprime 777-1,7712), and we can calculate the 
numbers f(m) if we know them for m = p r with p prime and r > 1. This simplification is 
a nice algebraic result that need no longer hold for the analogous problem applied to planar 
modules with iV-fold symmetry when N > 46, see |40) . 

In our present case, due to multiplicativity, a Dirichlet series ^(s) is an appropriate gen- 
erating function ^15^]. To calculate f(m) explicitly, we observe that, if m = p r is a prime 
power (p = 1 (mod 4), r > 1), only the two rotations 

( 

\U p 

lead to CSLs of index m, wherefore we have f(p r ) = 2 in this case. Now, we can directly 
calculate the entire generating function through its Euler product representation and obtain: 

Proposition 3.1. The Dirichlet series generating function of the number f(m) of CSLs of 
Z 2 of index m is 

•w = f^= n 

m=l psl(4) P=l(4) 
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A little later, we shall express ^(s) in terms of (^-functions. This generating function is not 
only a succinct way of representing the statistics of CSLs and coincidence indices, it is also 
a powerful tool for determining their asymptotic properties. For example, it can be used to 
show (through the determination of the right-most pole of <P(s) in the complex s-plane, which 
is at s = 1, and its residue) that the number of CSLs of Z 2 with index < N is asymptotically 
N/tt (and the corresponding number of coincidence rotations is asymptotically 4N/ir). The 
possible coincidence indices are precisely the numbers m with all prime factors = 1 (mod 4), 
and we then have 

(3.6) /(to) = 2 a , 

where a is the number of distinct prime divisors of to. Each CSL is itself a square lattice 
|4()[ 15]. with the index as the area of its fundamental domain. We shall come back to this in 
a more general context. 

Finally, the full group of coincidence isometries, OC(Z 2 ), is the semidirect product of the 
rotation part SOC(Z 2 ) (normal subgroup) with the cyclic group C 2 generated by complex 
conjugation (= reflection in the x-axis): 

(3.7) OC(Z 2 ) = SOC(Z 2 ) x C 2 . 

Here, conjugation of a rotation through an angle (p by complex conjugation results in the 
inverse rotation (through —(p). Let us give a justification of Eq. (|3.7|) . Since 0(2) = SO(2)xC2 
(semi-direct product) with the C 2 of Eq. (j3.7j) . any planar isometry T with det(T) = — 1 can 
uniquely be written as the product 

(3.8) T = R(<p) ■ T x 

of a rotation through ip with T x , the reflection in the x-axis (complex conjugation). But T x 
leaves the entire lattice Z 2 invariant, whence T is a coincidence isometry if and only if R(tp) 
is a coincidence rotation. 

The calculation of coincidence indices is also simple in this case. The coincidence index for 
the reflection T x is 1. For an arbitrary element of OC(Z 2 ), we either meet a rotation (where 
we know the result already) or use the factorization of Eq. ()3.8j) again. Then, the coincidence 
index is identical with that of its rotation part, so Eq. (|3.8j) is all that is needed. This solves 
the coincidence problem for the square lattice completely and we have 

Theorem 3.3. The group of coincidence isometries of the square lattice Z 2 is 

OC(Z 2 ) = 0(2, Q) ~ (C 4 x Z^) x C 2 • 

This group is fully characterized by Eqs. (|3,3[) . (|3.4|) and (|3.7j) . The coincidence index of a 
rotation is given by Eq. ()3.5|) , while that of a reflection (|3,8|) equals the index of its rotation 
part. The corresponding Dirichlet series generating function is ^(s) from Provosition \S.l\ □ 

3.2. A short digression on a hierarchy of problems. It is the intention of this paragraph 
to shed some more light on the coincidence problem and how it relates to similar questions. 
We shall explain it for the square lattice in an informal manner. 

To this end, let us start with the question of how many sublattices of Z 2 have index m - 
without any further restriction. Let us call this number a m . Clearly, a 1 = 1 (only Z 2 itself 
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is sublattice of index 1) and a 2 = 3 (counting two different rectangular sublattices and one 
square sublattice). In general, a mn = a m a n when m, n are coprime, and one can derive, either 
from the Appendix or from |51l Lemma 2 on p. 99], the general result 



^lH = lld\m d 



with generating function 
On 

(3.9) ~i w 



f{s) = Y.—s = as)- as -i) 

^— ' rn. s 



— i I 3 | 4 , 7 | 6 I 12 , 8 | 15 I 13 | 18 | 12 , 28 | 14 , 

t 2 s 3 s 4 s 5 s 6 s 7 8 s 9 s 10 s 11 s 12 s 13 s 

Here, £(s) = X)™=1 m ~ S * s Riemann's zeta function, compare [SO]. Prom this, it can be shown 
that the number of sublattices with index m grows, on average, linearly as mit 2 /6, compare 
|301 Thm. 324]. These results are, of course, affine in nature and apply to any lattice of rank 
2, and to any free Abelian group of rank 2 (counting the different subgroups of index m). 

Let us look for metric properties by asking how many of the sublattices of Z 2 of index m 
are actually square lattices (see (U El El for various generalizations of this question). This 
number can be obtained by counting the lattice points on circles of radius m (hence counting 
solutions of the Diophantine equation x 2 + y 2 = m) and afterwards dividing by 4 (the order 
of C4, the rotation part of the point symmetry group of Z 2 ). The result is given in Chapters 
16.9, 16.10 and 17.9 of [3U] and leads to the Dirichlet series generating function 

F (s) = CM = Y^~2r* II (t^) 2 - II TZp* 

(3.10) ^1(4) Vi P J P ^(A) L P 

= 1 + — + — + — + — + — + — + — + — + — + — + — + — -I 

2 s 4 s 5 s 8 s 9 s 10 s 13 s ~ 16 s ~ 17" ~ 18 s ~ 20 s ~ 25 s ^ 

where here and in what follows Ck( s ) * s the Dedekind zeta function of the quadratic field 
K = Q(i), compare |491 § 63, A. 14 on p. 251]. The average value of the coefficients of F(s) is 
constant, namely vr/4, which follows either from the asymptotic properties of the generating 
function near its right-most pole (at s = 1) or just from counting one quarter of the lattice 
points inside the circle of radius \/rn, which is mix /4 to leading order in m. 

In the last case, some of the square lattices fail to be primitive (such as 3Z 2 etc.), i.e., 
whenever the sublattice is an integer multiple of Z 2 or one of its primitve sublattices. If we 
exclude those, primes p = 3 (4) are impossible as divisors of the index m, and some solutions 
of p = 1 (4) also drop out (whenever the index is divisible by a square). Now, the generating 
function reads 

F( S ) = (i + 2-). TT i±£l = %4 

(3.11) psl(4) y SK ' 

= 1 + — 4- — -I — - — I — - — I — - — I — - — I — - — I — - — I — - — I — - — I — - — I- • • • 

L ^ 2 s t 5 s T 1Q 3 T 13 s -T 17 s T 25s ' 26 s ^ 29 s ^ 34 s T 37 s T 41 s ^ 

and the average number of primitive square sublattices of index m is given by 3/(2-7r). This 
can be determined by counting one quarter of the visible points .L in the circle of radius y^rn, 
which is (7rm • 6/7r 2 )/4. 
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Still, not all primitive square sublattices are CSLs - they only are if the index is odd. This 
finally results in the generating function of the coincidence problem described above, namely 

(3.12) #w = n ^ = ( i+2 " s )"|S 

p=l (4) 

where the average number of CSLs of index m is 1/tt. This equation expresses the Dirichlet 
series generating function in terms of zeta functions. Similar formulas will also appear in later 
examples. We hope that this short digression has put the problem in a broader perspective. 
Let us now climb up to higher dimension where the picture changes significantly because 0(d) 
is no longer Abelian for d > 2. 

3.3. d = 3: the three cubic lattices. In this paragraph, we shall use the notation -Tf,p,s 
for the F-type {fee), P-type, and Z?-type (bee) cubic lattices, respectively. In particular, 
Pp = Z 3 . Let us start with this and write the lattice as Z 3 = Zei © Ze2 © Ze^. We already 
know that OC(Z 3 ) = 0(3, Q) (other isometries might also lead to coincidences, but not to a 
CSL of full rank 3). The subgroup of rotations with index 1 is the rotation symmetry group 
of the cube of order 24, Aut+(Z 3 ) = = SO(3,Z), cf. 2 for details. 

As is well known, 0(3) = SO(3) x C2 is a direct product, where C2 = {±13} is the centre 
of 0(3). Consequently, we can restrict our attention to pure rotations, as reflections may 
be written as a product of a rotation R with — 1 3 (note that — Z 3 = Z 3 ). At this point, 
we introduce quaternions q, see [121122100], and Cayley's parametrization with 4 real 
numbers (k, A, fi,u) = q / 0, 

j / k 2 + A 2 - u 2 - v 2 -2kv + 2Xfi 2ku + 2\v 

(3.13) R(q) = —[ 2kv + 2\h k 2 — \ 2 + fi 2 — v 2 -2n\ + 2^v 

V -2npL + 2\v 2k\ + 2[iv k 2 - A 2 - [i 2 + v 

where |g| 2 = k 2 + A 2 + /x 2 + v 2 is the so-called reduced norm of the quaternion q. In particular, 
the multiplicative group of quaternions of norm 1, which form the unit sphere S 3 , provides the 
usual double cover of the rotation group SO(3,K), via the group homomorphism q \— > R(q) 
(since R(q) = R(—q)). 

As we are interested in rotation matrices with rational entries (i.e., in the subgroup 
SO(3,Q)), we consider Cayley's parametrization with 4 integers (k, A, /x,z/) = q / 0, and 
choose them coprime, i.e., gcd(/-t, A, /i, v) = 1. We call such quaternions primitive. This way, 
we parametrize the entire group SO(3, Q), and obtain each element of it exactly twice (again 
because of R(q) = R(—q)). With this approach, one obtains the following result. 

Proposition 3.2. Consider r p = Z 3 and let R = R(q) £ SO(3, Q) be parametrized by a 
primitive quaternion q. Rs coincidence index £(R) is the denominator of R as defined in 
Eq. (EU). It is the "odd part" of \q\ 2 , 

S(R) = den(R(q)) = \q\ 2 /2 e , 

where £ is the largest integer such that 2 e divides \q\ 2 . 
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Before we prove the statement, let us remark that the necessity for the division by 2 e 
in Proposition 13.21 stems from the fact that some primitive quaternions (such as (1,1,0,0) 
or (1, 1, 1, 1)) have norms that are not coprime with the matrix entries of R{q). A closer 
inspection, using results of [32], shows that only the exponents t G {0, 1, 2} are possible. 

Proof: Consider R(q) with ^ q = (k, A, fi, v) primitive and a = \q\ 2 . Define the 4 vectors 

v = I H I , »i = I v I , v 2 = 

which all have integer preimages under R = R(q), and the 4 matrices 

B o = ( v i> v 2i v 3) > det( J B ) = kct ; B x = (u , t> 2 , « 3 ) , det(B 1 )=Xa; 

B 2 = ( v i, v o, v 3) i det( J B 2 ) = //ct; 5 3 = (t^, « 2 , « ) > det(5 3 ) = z^r . 

Now, each i?j can be read as a basis matrix of a lattice that is a sublattice of both Z 3 and 
i?Z 3 , hence it is also a sublattice of the CSL 1? n i?Z 3 . Consequently, the coincidence index 
A7(i?) must divide each of the determinants det(J3 3 -). As q was primitive, S(R) must therefore 
divide a. Since S(R) trivially also divides the third power of den(i?) = a/2 e (which is odd), 
we have E(R)\dea(R). 

To establish the claim, it is now sufficient to show that also den(R)\U(R). Observe that 
[I? : I? n Rl?\ = [RI? : I? n Rl?\ due to det(i?) = 1, as in the proof of Lemma l2~H Since, 
by definition of the denominator, 

gcd{den(i?) J R ij - | 1 < i,j < 3} = 1 , 

the lattice i?Z 3 contains a vector of the form a/den(i?) with a primitive vector a £ Z 3 , i.e., 
a = (a 1; a 2 ,a 3 )* with d; 6 Z and and gcd(a l5 a 2 , a 3 ) = 1. This implies that the number of 
cosets of Z 3 n i?Z 3 in i?Z 3 must be a multiple of den^), so that den(R)\U(R). □ 

In particular, this reproduces the well-known result [21J EHj that E(R) = E z a(R) runs 
precisely through all odd integers, i.e., A7 Z 3(SO(3, Q)) = 2N + 1 = {1, 3, 5, 7, . . .}. 

Cayley's parametrization has the nice property that v Q = (A, fj,, v) 1 gives the (generic) 
rotation axis of R(q) = R(k, A, fi, v), 

(3.14) R(k,\,h,u)v = v , 

while the rotation angle follows from tr(R) = 1 + 2cos((£), so 



(3.15) cos(99) 



k 2 + A 2 + /U 2 + 



2 



One can easily construct all solutions for small indices explicitly, while the case k = gives 
all coincidence rotations through tt as described by Luck [31]. If 24/ (m) is the number of 
coincidence rotations of index m, the arithmetic function f(m) once again counts the different 
CSLs of index m. This function is multiplicative as a consequence of the fact that integer 
quaternions 1 have unique left (and right) factorization up to units 2 , see JSH] for details. 

lr They constitute the Hurwitz ring JT of quaternions of the form \{a, b, c, d) with a, b, c, d all even or all odd. 
2 They are the 24 quaternions ±(1,0,0,0) (and permutations) and |(±1, ±1, ±1, ±1). Together, they form 
a group that is the double cover of the symmetry group (rotations only) of the regular tetrahedron, see 0. 
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Its calculation amounts to counting the representations of an integer as a sum of 4 squares, 
compare 32, Ch. 11] or |3()| Ch. XX], and to observe the relation between \q\ 2 and den(R(q)). 
This results in 

/(l) = 1 and /(2n) = 0, 

f(p r ) = (P+l)p r 1 1 f° r °dd primes and r > 1, and 
f(mn) = f(m)f(n) , for m,n coprime (multiplicativity of /), 
see also [261 158| . This can be summarized as follows. 

Proposition 3.3. The Dirichlet series generating function for the number f{m) of CSLs of 
Z 3 of index m reads 

fOr-.y* /(m) - TT 1+p ~ s = 1 - 2l " c(«)C(«-i) 

(3.16) ^ ™° ~ l^-P 1 '* 1 + 2 ~ S C(2«) 

- 1 I 4 I 6 I 8 I 12 I 12 I 14 I 24 I 18 I 20 I 32 I 24 I 30 I 
* r 3 s " r 5 s " r 7 s 9 s " r 11 s " r 13 s " r 15 s 17 s " r 19 s ^ 21 3 23 s 25 s 

With this generating function, one can again determine the asymptotic behaviour of f(m). 
The result is that the number of CSLs of index < N is asymptotically given by 3N 2 /ir 2 , and 
the number of coincidence rotations with index < N by 72N 2 /it 2 . 

To expand on the systematics of our generating functions, let us remark that Eq. ()3.16j) 
can be rewritten as 

(3.17) 0(s) = • 

where Ch( s ) = (1 — 2 1_2s ) C(2s) C(2s — 1) is the zeta function of the Hurwitz ring J of integer 
quaternions, compare |49l § 63, A. 15 on p. 252]. It is the generating function for the number 
of non-zero right ideals of the ring J, which is a maximal order of the standard quaternion 
algebra over Q. 

This is not quite the end of the story. Coincidence rotations of a given index m can be 
collected into equivalence classes of rotations related by the action of the point symmetry 
O or, more generally, of 0(3, M). This requires a double coset analysis that is described in 
|261 158j . It turns out that inequivalent CSLs of 1? with the same index occur for the first 
time at E = 13. 

Also, describing the fine structure of a coincidence rotation requires an analysis of the 
lattice planes perpendicular to the rotation axis. For example, the (unique) equivalence class 
for S = 3 can be represented by q = (0, 1, 1, 1), i.e., a rotation through tt around (1, 1, 1)*. 
Here, three layers are stacked periodically (with period v / 3), with perfect coincidence in one 
(called A, which is therefore an ideal grain boundary), but none in the other two (called B, C, 
which are interchanged in type by the rotation). So, this coincidence rotation can be used to 
create a twin that shows up as a change in the stacking sequence, see Figure ^ 

Let us take a short look at the other cubic lattices, Tp and Tb (from now on, we use indices 
F,P,B to distinguish the three cubic lattices). They possess the same group of coincidence 
isometries as Jp = Z 3 , 

(3.18) oc(i» = oc(r P ) = oc(r B ) = 0(3, q), 
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Figure 1. Construction of a twin based on R(q) with q = (0, 1, 1, 1). 
but even more, one has 

Proposition 3.4. In E 3 , i/ie coincidence index of an orthogonal matrix is the same for all 
three cubic lattices. 



Proof: We can apply Lemma 12.31 From the inclusion 
(3.19) 



2 2 

r F c r P c r B , 



we know that Up \ 2Up and Up \ 2Up. But Up = Up from Theorem 12.21 because P^ 
in this setting. Since Up is always odd, we can either have Up = Up or Up = 2Up. 

However, as Tp and Pp have the same point symmetry, Up = 1 implies Up = 1, and Tp 
cannot possess coincidence isometries of index 2. From the multiplicativity of the number of 
CSLs of index m, see [581 HUj . it then follows that no CSL of r B can have even index, hence 
Up = Up = Up. D 

One can see the second step also without reference to the multiplicativity of the CSL 
counting function. Observe that Tp = fpU (v + 7p), with v = 1, 1)*. These two cosets 
have disjoint shells, and one can check that coincidence rotations produce coinciding points 
in both shells with equal density, which then implies that Up = Up. 

Nevertheless, there are specific differences in the coincidence structure of the three cubic 
lattices. They show up in different layer arrangements [2E], but details go beyond the scope 
of this article. 



3.4. d = 4: the two hypercubic lattices. In four dimensions, there are only two hyper- 
cubic lattices (P-type and centred, represented by Z 4 and D4), and they confront us with a 
different situation that occurs in no other dimension: the centred lattice, D4, has a larger 
symmetry than its partner of P-type, Z 4 . The dual lattice of D4, the weight lattice D%, is 
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equivalent to D± (through a rotation followed by a homothety) and need not be considered 
separately, compare [T%] . 

Let us start with the description of the coincidence structure of -D4. First, we observe that 
quaternions give us again a helpful parametrization of rotations |22l I33| . where we now need 
a pair {q l ,q 2 ). The corresponding matrix is defined through 

M(q 1 ,q 2 )x t := {q-^xq^f , 

where x t stands for the transpose of the quaternion x (and is thus a column vector). With 
q 1 = (k,£,m,n) and q 2 = (a,b,c,d), one finds the 4x4-matrix 



M( qi ,q 2 ) 



/ ak+b£+cm+dn —a£+bk+cn — dm ~am—bn+ck+d£ —an+bm—c£+dk\ 

a£—bk + cn — dm ak+b£—cm—dn —an+bm+c£—dk am+bn + ck + d£ 

am—bn—ck+di an+bm + c£+dk ak—b£+cm—dn —a£—bk+cn+dm 

\an+bm — c£~dk —am+bn—ck+dl a£+bk + cn+dm ak — b£ — cm+dn/ 



which has 

dek(M(q u q 2 )) = {k 2 + I 2 + m 2 + n 2 ) 2 ■ (a 2 + b 2 + c 2 + d 2 ) 2 = |Qi| 4 |g 2 | 4 

and satisfies the orthogonality relation 

MM 1 = ^det(M) ■ 1 4 . 

Consequently, |qr 1 | = \q 2 \ = 1 results in a 4d rotation matrix, and the group homomorphism 
M: S 3 x§ 3 — > SO(4) is onto. It provides a twofold cover of the rotation group [S3]) with 
M(q 1 ,q 2 )=M(-q 1 ,-q 2 ). 

From here, we can find a parametrization of SO(4, Q) if we start from two non-zero primitive 
integer quaternions q^q 2 (i.e., each quaternion has the form (k,\,(j,,v) with k,\,(j,,v £ Z 
and gcd(«;, A, fi, v) = 1). Then, if we consider the matrix 

(3.20) R(q l ,q 2 ) ■= Mf^-,^ 2 -) = -J—M( qi ,q 2 ), 



ml ml/ mQzl 

we see that it is an element of SO(4, Q) if and only if |QiQ2l 2 ^ s a square in N, in which 
case we call the pair of integral quaternions admissible. But with all admissible pairs of 
primitive quaternions, we actually exhaust SO(4, Q), and obtain each element twice (due to 

RfavQi) = R(-Qi,-Q2))- 

From Theorem 13.11 we already know that (S)OC(Z?4) = (S)0(4, Q). As in the previous 

cases, it is sufficient to treat rotations, since reflections can be written as a product of a 

rotation with a special reflection that leaves D4 invariant - in complete analogy to the situation 

in the square lattice (0(4) = SO(4) x C2 is a semi-direct product). So, we need to know the 

coincidence index of an arbitrary rotation R G SO(4, Q). This is not just the denominator of 

R, but given by the following result, see [591 114j for a proof. 

Proposition 3.5. Let (q 1 ,q 2 ) ^ e an admissible pair of primitive integral quaternions, and 
let S{q) denote the index defined above in Eq. (|3,2[) . Then, the matrix R(q 1 ,q 2 ) £ SO(4, Q) 
has coincidence index 



(3.21) 
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Here, 1cm denotes the least common multiple, and the subscript F refers to the (face-) 
centred lattice D4. With this formula, it is now a combinatorial problem to determine the 
number of coincidence rotations of a given index m and, dividing by 576 (the order of the 
rotation symmetry group of -D4), also the number f F {m) of different CSLs of -D4 of index 
m. As follows explicitly from Eq. Q3.21JI and also from the unique factorization of integral 
quaternions, f F (m) is again a multiplicative function, so that we only need to calculate it for 
m a prime power. 

Starting from Eq. ()3.21j) and counting the possibilities to contribute to U F (p r ), it is not 
difficult to derive (for r > 1) the explicit expression 



[§] 

-2b 



(3.22) f F (pl = f(p r ){f(p r ) + 2J2f(P r 

e=i 

with the f(m) of the 3d cubic case in Eq. (|3.16|) . An empty sum is to be understood as 0, 
and [.] denotes Gaufi' brackets. The result is (see (55] for a proof) 

/f(1) = 1, 
f F (mn) = f F (m)f F (n), if m, n coprime , 
f F (2m) = 0, and 

f F (p r ) = — p r ~ 1 (p r+1 + p r ~ 1 — 2) , for odd primes and r > 1 . 
p — 1 

This fixes the Dirichlet series generating function and one obtains 

Proposition 3.6. The Dirichlet series generating function for the numbers f F (m) of CSLs 
of index m in the root lattice D 4 reads 

g n _ y /fH (i +p -*)(i +p i- 8 ) 

(3.23) F{> h ™* li ( l-pi-*)(l-p2- s) 

— 1 I 16 I 36 I 64 1 168 , 144 , 196 , 576 , 324 . 400 , 1024 , 576 , 

— 1 T js T js T js T gs t ip T 13 s ~T 15 s ~T yja "t" lgs ~T 2 1>> ~l~ 23 s ^ ' ' ' 

Again, the possible indices are precisely all odd integers, Z , £> 4 (SO(4, Q)) = 2No + 1. A 
comparison with the cubic case in 3-space reveals the remarkable identity 

(3.24) f (b) = *(a)#(a-l) 

where <&(s) is the generating function of Eq. (|3.16j) . This makes is rather easy to calculate 
the asymptotic behaviour from that in 3 dimensions. The right-most pole of @ F (s) is at 
s = 3, so one obtains that the number of CSLs of index < iV grows asymptotically as 
^£(3)iV 3 ~ 0.26257 iV 3 (note that ((3) is known to be irrational, but its value is only known 
numerically). For the number of coincidence rotations with index < N, one has to multiply 
by 576. Another consequence of Eq. (|3.24j) is the formula 

(3.25) f F (m) = 5>-/(d)./(m/d) 

d | m 

which follows from the convolution theorem of Dirichlet series and allows for a simple and 
efficient calculation of the numbers f F (m). 



THE COINCIDENCE PROBLEM FOR d < 4 



17 



Having described the root lattice -D4, the (face-) centred cubic lattice in four dimensions, 
we now turn to the slightly more complicated case of the P-type cubic lattice, Z 4 . Prom the 
inclusion 

(3.26) D A c Z 4 c D\ 

and the result that D 4 and D\ have the same index formula, see Theorem 12 .21 we get Up \ 2Up 
and Up I 2Up. But Up is always odd, so either Up = Up or Up = 2Up. Here, in contrast to 
the situation in 3-space, both possibilities arise. The point symmetry group of D4 is larger 
than that of Z 4 , with 

(3.27) [Aut(D 4 ) : Aut(Z 4 )] = 3. 

Consequently, one third of the elements of Aut(-D 4 ) are symmetries of Z 4 while the others 
result in coincidence isometries of Z 4 . They turn out to have index 2. In going from here to 
the number of CSLs of index m, it is clear that we have /p(l) = 1, /p(2) = 2 and f p (2 r ) = 
for r > 1. The multiplicativity (which needs to be proved, e.g., similarly to the arguments 
given in |59| ) of fp{m) then gives the general answer 

/p(m) = fp(m) , for m odd, 
f P {m) = 2f F (m/2) , for m = 2 (4), and 
/ P (4m) = 0. 

This can now easily be summarized as follows. 

Proposition 3.7. The Dirichlet series generating function for the number fp(m) of CSLs 
of index m in Z 4 reads 



(3.28) i 2 ( 1_p s )( l ~P 2 s ) 



— ~r 2 s 3 s 5 s 6 s 7 s 9 s * r 10 s ~ l ~ 11 s 13 s ~ l ~ 14 s 15 s 17 s ~ l ~ ' ' ' 

Note that, in comparison to the -D4 case, the number of CSLs grows faster by a factor of 
5/4, while the number of coincidence rotations grows slower by a factor of 5/12, due to the 
smaller point symmetry group of Z 4 . 

This calculation was actually possible without giving the corresponding index formula first, 
by using the multiplicativity of f p (m). Since SOC(Z 4 ) = SOC(L>4) = SO(4, Q), there must 
be a different index formula for Z , and indeed one obtains, for an admissible pair of primitive 
integral quaternions (see [25] for a proof) 

U P (q 1: q 2 ) =lcm {U{q x ), U(q 2 ),den(R(q 1 , q 2 ))} 
=\cm{U F (q 1 ,q 2 ),den(R(q 1 ,q 2 ))} . 



At this point, we close our description of the coincidence structure of lattices and turn to the 
perhaps more interesting case of quasicrystals. 



18 



MICHAEL BAAKE 



4. Coincidence isometries for modules 



So far, we have described the case of lattices and crystals. For the treatment of quasicrys- 
tals, we have to extend our concepts to Z-modules, embedded in Euclidean space, which are 
not necessarily discrete any more. 

Definition 4.1. A subset Ai of E d is called a ^-module, of rank r and dimension d, when it 
is the Z-span of r vectors a 1 , . . . ,a r (the basis of the module) that are linearly independent 
over Z, but span E rf over R. 

Clearly, we must have r > d, and a module with r = d is a lattice. As a group, a module of 
rank r is isomorphic to the free Abelian group of rank r - so we may consider such modules 
as special geometric realizations of free Abelian groups. 

At this point, the concepts of submodule, commensurateness, and that of coincidence isom- 
etry and index, are defined in exact analogy to Section [2j so there is no need to repeat them 
here. We only note that the set of coinciding points is now a module rather than a lattice, 
wherefore we call it a coincidence site module, or CSM for short. A bit more care is needed 
for the definition of a dual module. Let us first observe 

Lemma 4.1. For every module M. C E d of rank r and dimension d < r, there is a lattice 
r C E r such that M. is the one-to-one projection of T into K d . 

Proof: Out of the basis {a 1 , . . . , a r } of A4, we can pick n R-linearly independent vectors 
which span E rf over R, say a 1; . . . , w.l.o.g. They span a lattice in K d . The statement is now 
obvious as we can add one new dimension for each basis vector remaining. This can always 
be done in such a way that the projection is one-to-one on r. □ 

It is an obvious idea to try to define a dual object for a module through a lift to a lattice 
r because then T* is well-defined and can be projected down again. Unfortunately, this is 
neither unique nor satisfactory, as it can happen that the object defined this way is a module 
of smaller rank than the original one. If, however, there is some additional structure (e.g., 
irreducible symmetry), such a lift can be made essentially unique and the dual object is well- 
defined, see |381 139j for relevant examples in our present context. If this situation applies, it 
is again true that a module and its dual module share coincidence group and index formula. 

In what follows, we concentrate on examples that are connected with the golden ratio, 
t = (1 + y/5)/2. The modules that will appear are invariant under multiplication by r or a 
power thereof. So, it might be instructive to see how the construction of a dual module works 
here. The simplest example is the ring of golden integers 



which is the ring of algebraic integers in the quadratic field Q(r), but can also be seen as 
a Z-module M of rank 2 and dimension 1 in E. The special structure that helps in this 
case is the existence of an automorphism ' of Q(r), called algebraic conjugation, which maps 
t to its algebraic conjugate r 1 = — 1/r = 1 — r (hence the name) and thus a = a + br to 
a' = a + br' = (a + b) — br. Now, the set 





(4.2) 
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is a lattice in 2-space E 2 , from which one obtains Z[r] by projection into E. In E 2 , r has a 
well-defined dual lattice, r*, and neither r nor r* has a lattice direction parallel to E. But 
then, we define the dual module M* by the projection of r* into E (which, in this case, gives 
M* = M/V5 = Z[t]/VE), and this object is unique in the sense that any other embedding 
of M into E 2 which maps algebraic conjugation to a lattice automorphism will result in the 
same dual object Ai*. 

Though examples in higher dimensions are more complicated (and will require more than 
just one automorphism), the basic idea is similar. In particular, it applies to all quasicrystals 
of interest, compare jHS]. It turns out that root lattices |18j prove extremely handy here UJ. 
Let us now illustrate the coincidence problem for noncrystallographic patterns by a seris of 
examples related to fivefold symmetry. For further material on planar structures with iV-fold 
symmetry, we refer to I4f)[ |2j . 

5. Modules and quasicrystals: the root systems H 2 , H 3 and H4 

Among the many possible quasiperiodic tilings, those attached to fivefold symmetry are 
of particular interest, especially in view of their application in solid state physics. Let us 
therefore start from the exceptional Coxeter groups H2 (usually called ^(S)), #3 and H$ |31| 
shown in Fig. [2j They are the symmetry groups of certain regular polytopes |19j . namely 
of the regular decagon ({10} in Schlafli's notation), the icosahedron {3,5} (or dodecahedron 
{5, 3}) and the regular 600-cell {3, 3, 5} (or the regular 120-cell {5, 3, 3}), and are of order 20, 
120 and 14400, respectively. 

C^-O O 0^<3 O O 0^-<D 

Figure 2. The non-crystallographic Coxeter groups of type H. 

The corresponding root systems (up to normalization) are given by the vectors which point 
to the 10 vertices of a regular decagon (H2), to the 30 vertices of the icosidodecahedron 
{ I } (i/3), and to the 120 vertices of the regular 600-cell in 4-space (H4), see (SHI or j3J for 
details. The Z-spans of these root systems define ff n -symmetric modules, which can be seen 
as projections from the root lattices A4, Dq and Eg. In particular, they all have well-defined 
duals, though we do not expand on this question here - it is discussed in detail in [3*5] . 

Having set the scene, we can now describe the coincidence structure of these modules (and 
some closely related ones). We shall also briefly discuss the connection to quasiperiodic tilings. 

5.1. d = 2: Coincidence rotations for tenfold symmetry. Let us consider a 2d qua- 
sicrystal with tenfold symmetry, the Tubingen triangle tiling 7 of Fig. |3J say. As mentioned 
earlier, the coincidence problem splits into two parts: first, the coincidence problem for the 
underlying Z-module M\q (which is the limit translation module ^5] of the tiling) and sec- 
ond, the correction, due to the acceptance domain, of the coincidence indices obtained in this 
way [HI]. Here, we discuss in detail only the first part. 
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Figure 3. Cartwheel version of the Tubingen triangle tiling. 

In the complex plane, the tenfold module of rank 4 over Z |35j can be written as the direct 
sum 

(5.i) Mio = ii e z^ e z^ 2 e z^ 3 , 

with £ = e 2,U//5 a fifth root of 1. As such, .Mio is the ring of algebraic integers in the cyclotomic 
field K§ = Q(£) (a field extension of Q of degree 4). Again, prime factorization is unique up 
to units 55 . One can now go through essentially the same argument as in the case of the 
square lattice. The primes in Z[£] are slightly more complicated, but only those dividing a 
rational prime p = 1 (5) can enter the factorization of e l(f = a/ (3 with a, j3 € Z[£] coprime. 
The result is 

Proposition 5.1. Every coincidence rotation (written as e llp £ K§ = Q(£)) of the standard 
tenfold symmetric module M\q can be factorized as 



(1) (2) 
(1) \ n v / (2)\ n P 



i / I n (2) 



(5.2) e*= e J] (3 

psl (5) 

where e is a Wth root of 1 and thus a unit in Z[£], and on/y finitely many of the exponents 
rip 2 ^ are different from 0. 

This factorization is slightly more complicated than that in the case of the square lattice, 
as we have two independent generators for each basic index p = 1 (5). This originates in the 
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fact that these (rational) primes are the product of 4 cyclotomic primes in Z[£] which form 
two independent pairs of complex conjugates. 
The index of a rotation R is 



(|4 13 |+|4 2) I) 



(5.3) U(R) = J] p 

p=l (5) 

and the group of coincidence rotations has the form 

(5.4) SOC(TWio) ~ Cio x . 

Thus, in spite of the more complicated factorization, the structure of the coincidence group 
remains simple (as it does for other planar symmetries, compare |4()1 [2]). Let us give two 
examples before we continue. U = 11 and E = 31 are the smallest non-trivial indices, with 
two generators each. With £ = e 2m ^ , they can be written as follows: 

2 + £ 2 + £ 2 2 - £ 2 - £ 2 

(5.5 £ = 11: -4, -^o? £ = 31: 1, V 

2+£' 2 +f 2-£' 2 -f 

Coincidence reflections can be described as products of rotations with the reflection in the 
x-axis, exactly as in the case of the square lattice, so we need not repeat that argument here. 

If 10/(m) denotes the number of coincidence rotations of index m, the multiplicative func- 
tion f(m) actually counts the number of different CSMs of index to. We then obtain, with 

QW =0( v / 5) = Q(£)nfc 



Proposition 5.2. The Dirichlet series generating function for the coincidence problem of the 
tenfold symmetric module of rank 4 in the plane reads 



/(to) -pr fl+p~ s \ 2 _ 1 C 



(5.6) £i ^ Ji( 5 ) Vl " PV 1 + 5 ' S C ®(r)( 2s 



^ = e^ = n , 

p=l(5) V 

4,4,4,4, 1,8, i ,4, i 



Q(0 [ 



11 s 31 s T 41 s T 61 s ^ 71 s ^ 101 s T 121 s ^ 131 s T 151 s ^ 181 s T 

Here, Cq(^)(' s ) is the Dedekind zeta function of the cyclotomic field Q(£), while Cq( t ) i s 
the zeta function of the maximal real subfield, Q(£ + £) = Q(t), see also Eq. (|5.1Up below. 

All CSMs are scaled versions of A^io and the number of CSMs of index < iV is asymptoti- 
cally 5 log(r)iV/7r 2 (while the number of coincidence rotations with index < iV is 10 times as 
large) . 

How do these results apply to the coincidence problem of the tenfold symmetric triangle 
tiling? The latter can be obtained through projection from the root lattice A4 with a regular 
decagon as its window [7j . A coincidence in the set of vertex points occurs if and only if there 
is a coincidence in the module Mio such that the image point in internal space lies both in 
the original window and in an appropriately rotated window. A consequence of this is that 
the coincidence group of the tiling is still SOC(A^io), but also that the £-factor or "index" 
of each group element is normally smaller than its index in Aiio by a correction factor close 
to 1 (depending on the group element). 

This is called the window correction factor, and has to be calculated for each tiling sepa- 
rately, as a function of the rotation angle. It also explains why the set of coinciding points 
forms a tiling of slightly different type from the original one, a small proportion of the points 
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of the original tiling being missing from it. In fact the term "index" for the reciprocal of the 
fraction of coinciding points is no longer appropriate in this setting, as it neither has a purely 
algebraic interpretation, nor is an integer any more. Details of this and the determination of 
the rotation angle in internal space by means of algebraic conjugation are given in ^Uj. In 
many examples, in particular when the windows are regular polytopes, the maximal error is 
so small (at most of the order of a few percent) that one can safely ignore it and work with 
the module index instead. 

5.2. d = 3: the icosahedral modules of rank 6. Icosahedral quasicrystals are of partic- 
ular interest, and one would like to know their coincidence structure in detail (541 141j . We 
restrict our discussion to the investigation of the 3 different 3d icosahedral modules of rank 
6 over Z |47j and again omit the determination of the window correction. We shall call the 
modules Mp, Mp, Mf for B-, P- and F-type, respectively 3 . They are spanned by the 
orthonormal basis ei,e2, £3 with coefficients on E Z[r], r = (1 + v / 5)/2, as follows: 

Mb = {Ef=i «i e i I r2 «i + T "2 + «3 = (2)} 
(5.7) M P = { x E M B I "i + a 2 + a 3 = or r (2)} 

Mf = { x G Mb I "i + a 2 + «3 = (2) } . 

The use of an orthonormal basis may be a bit surprising at first sight, but it will prove useful 
in a moment. It is possible in this simple form because e 1 ,e 2 ,e 3 are chosen parallel to 3 
mutually orthogonal twofold axes of the icosahedron. In this setting, Mf is the Z-span of 
the root system of type H3, and hence a Z-module of rank 6 and dimension 3. It actually also 
is a Z[r]-module (of rank 3), which is also true of Mb, but not of Mp. This relates to the 
fact that M f and M p are invariant under multiplication by r, while Mp is only invariant 
under multiplication by r 3 . 

To describe the coincidence rotations (in the orthogonal basis), Cayley's parametrization 
can again be used. Our first assertion is that the coincidence group is the same for all three 
modules: 

Proposition 5.3. OC(M B ) = OC(M P ) = OC(M F ) = 0(3, Q(r)) . 

Proof: Observe the relations 2Z[t] 3 C M f C M p C M f C Z[r] 3 and Z[t] 3 = Z 3 tZ 3 . 
These modules all possess the same OC-group, and this obviously is 0(3, Q(r)) by the same 
type of argument we have used previously in the discussion of the (hyper-)cubic lattices. □ 

The unit quaternions (1, 0, 0, 0), |(1, 1, 1, 1), |(r, 1, — 1/r, 0) together with all even permu- 
tations and arbitrary sign flips form a group Y of order 120 which is the usual double cover 
(HUGH of the icosahedral group Y = {R G SO(3,Q(r)) | U(R) = 1}. The icosian ring I, see 
|36j for details, consists of all integral linear combinations of elements in Y and is a maximal 
order with unique (left- or right-) factorization in the quaternion algebra over the field Q(t). 
One finds the relation SO(3,Q(r)) = {R(q) |0 ^ q G I}, and our second assertion is the 

This terminology originates from the fact that these modules can be obtained as projections of the three 
types of hypercubic lattices in 6-space, Dg, Z 6 , and D 6 . 
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index formula for a coincidence rotation Rq £ SO(3,Q(r)), again for all three modules: 
(5.8) U(R ) = gcd{iV(|q| 2 ) | g€l, R(q) = R } , 

where the argument \q\ 2 on the right hand side is always a number in Z[r] and its norm is 
defined by N(m + nr) = m 2 + ran — n 2 . We use the convention that the gcd is always a 
positive number. The indices E run through all positive integers of the form m 2 + ran — n 2 
with integral m and n. These are the numbers all of whose prime factors congruent to 2 or 
3 (mod 5) occur with even exponent only. (They can also be characterized as the positive 
numbers of the form 5x 2 — y 2 with integral x and y, as used in For U < 100, one finds 

the list of numbers 

1, 4, 5, 9, 11, 16, 19, 20, 25, 29, 31, 36, 41, 44, 45, 49, 55, 59, 61, 64, 71, 76, 79, 80, 81, 89, 95, 99, 100, 

which covers the cases known from |54j . 

If 60/(m) is the number of coincidence rotations of index m, f(m) is the number of different 
CSMs of index m. Since the icosian ring is a maximal order with unique (left- or right-) 
factorization |461 152j . f(m) is again a multiplicative function, i.e., /(I) = 1 and f(ran) = 
f(m)f(n) whenever m,n are coprime. Furthermore, with r > 1, one finds 

/(5 r ) = 6-5 r -\ 

Then, if p = ±2 (5), 

/( p 2r-i) = and f{p 2 r) = (p 2 + 1)p 2(r-l)_ 

Finally, if p = ±1 (5), 

f(pl = (p+l){{r + l)p r - l + {r-l)tf- 2 ). 
This fully determines the generating function of f{m). 

Proposition 5.4. The Dirichlet series generating function for the number of CSMs of an 
icosahedral module from Eq. (|5.7j) is given by 

, n _ y /M _ 1 + -tt l+P~ 2s i-r f l+p-' V 

2^ m s 1-51-* 11 l_p2(l-«) 11 \l_pl--i 

m=l p =±2 (5) y p=±l (5) V ^ / 

- i I 5 | 6 | 10 | 24 , 20 i 40 , 30 I 30 I 60 , 64 | 50 I 

4 s 5 s 9 s 11 s 16 s 19 s ~r 20 s " r 25 s ' 29 s ' 31 s 36 s " r ' ' ' 

The number of CSMs of index < iV is asymptotically 45\/5 log(r)A^ 2 /27r 4 (while the number 
of coincidence rotations with index < N is 60 times as large). 
The function 'P(s) can be expressed in terms of zeta functions as 

(5.9, ,(,) = ™^ - 77TT 

with the quadratic field L := Q(r) = Q(\/5), Ci(s) = Cl(2s)Cl( 2s ~ the C-function of the 
icosian ring, and 

< 5 ' 10 ' &W-T^ II rr^ II (T^. 

p=±2 (5) ^ p=±l (5) V ^ ; 
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A more complete description of the icosahedral case, which requires a certain level of mathe- 
matical machinery and various results from algebraic number theory, will be given in 

5.3. Short digression on a related cubic module of rank 6. Crystals and quasicrystals 
are specific examples of ordered phases with longe-range order, but there are many other 
ones. Incommensurate structures are also widely studied in the literature, and they are also 
connected to modules rather than lattices. It is an obvious question what happens to such 
modules if they still have cubic symmetry (and hence dimension 3), but rank 6. There are 6 
types of such modules, called B + B, P + P, F + F, B + P, B + F, and P + F in a suggestive 
notation. Let us consider the case P + P in a bit more detail. It is clear that Z 3 + aZ 3 is an 
example of it if we demand a Q (otherwise, the rank would not be 6). 

So, let us assume that a is irrational. Then, Z 3 n aZ 3 = {0} and we can write the module 
as Z 3 © aZ 3 . Generically, the OC-group will be that of Z 3 itself (with the index squared), 
because, if a is not algebraic, there is no rotation which brings a point of Z 3 into coincidence 
with one of aZ 3 . This changes quite a bit if a is algebraic. A case of particular interest in 
our context is that of a = r, where we get the module 

M c = Z 3 © rZ 3 = Z[r] 3 . 

This module has the same OC-group as the three icosahedral modules above, OC{Mc) = 
0(3, Q(r)), but a different index formula. In fact, in complete analogy to the cubic lattices, 
one can prove that 

(5.11) S(R) = \N(den(R))\ . 

Here, den(i?) is defined w.r.t. Z[r] and hence a number of the form m + nr, N is the norm in 
Z[t] as used above, and the absolute value is needed because the denominator is only defined 
up to units in Z[r], which are the numbers ±r r with norm (— l) r . 

It is clear from Eq. (|5.11|) that the set of indices is the same as in the icosahedral case, i.e., 
all positive integers which are representable by the integral quadratic form m? + ran — n 2 . 
More surprising is the result that the generating function is very similar: differences only 
occur for indices which are divisible by 4. Explicitly: 

M*) = 1 + 4l ~* • Cl{s)Cl{s ~ 1} 

(5.12) C{ ' 1 + 4- ( L (2s) 

— i I 8 I 6 , 10 , 24 , 32 , 40 , 48 , 30 I 60 , 64 , 80 , 

— 4 s 5 s 9 s 11 s 16 s 19 s ~r 20 s 25 s " r 29 s 31 s 36 s " r ' ' ' 

Due to the prefactor, the number of CSMs of index < grows faster than that of the 
icosahedral case by a factor of 10/9, while, due to the different point symmetry groups, the 
number of coincidence rotations grows slower by a factor of 4/9. 

5.4. d = 4: the icosian ring as ^-symmetric module. Our final example is the icosian 
ring itself, compare |36j . viewed as the module obtained as the Z-span of the root system of 
type Hi. It is the limit translation module of a highly symmetric 4d quasicrystal studied in 
[2*3] . and can also be obtained by projection of the root lattice E% to a 4d subspace that is 
invariant under the action of the symmetry group of the regular 600-cell (which, in turn, is 
isomorphic to the Coxeter group H4), see [HlEZl for more. 
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The solution of this case proceeds in close analogy to that of -D4, details will appear in 
|14j . The coincidence group obviously is OC(I) = 0(4, Q(t)), and in order to obtain a 
parametrization of its rotation matrices, we are again using admissible pairs of quaternions 
g 1 ,g 2 £ I- This now means that | ^1^2 1 2 mus t be a square in Z[r]. This gives some slight 
extra complication from the arithmetic of Z[r], compare as rational primes p = ±1 (5) 
split into a product of two Z[r]-primes which are algebraic conjugates of one another, but 



not associates. Consequently, Eq. (j3.22j) has a counterpart in the present case that has to 
be modified for such primes by an extra factor 1/2 on the right hand side, while it is true in 
its unaltered form for all other primes. The result is a multiplicative function fJjn) which 
counts the CSMs of I of index in. It is specified by /j(l) = 1, and, for r > 1, by 

/n(5 r ) = |^- 1 (5 r+1 + 5'- 1 -2); 

if p = ±2 (5), one has 



f 1 (p 2r - 1 ) = and Mp 2r ) = ^P 2(r - 1} (P 2{r+1} +P 2(r - 1} -^)- 

p z — 1 ' 

Finally, for p = ±1 (5), one obtains the rather lengthy expression 

fl&) = {P [ p - V lf (4p 2 (2(p 2 + l)+r(p 2 -l))+/(p 2 + l)(r(p 4 -l)+p 4 -4p 3 -2p 2 -4p+l)). 
Although this looks a bit nasty, the corresponding Dirichlet series has the nice form 

<f>j( s ) = #( s )<p( s - 1) - ^i^W 



(5.13) C L (2s)C L (2s-2) 

1 i1j.1j_Mx18i!«j.8!»i?1i?60x 1800 , 2048 , 2500 , 
-L-|- 43 -|- 5s -|- gs T lp + Igs t 19 s "T 20" 25 s " r 29 s " r 31 s ~r 36 s " r ' ' " 

which resembles the situation of the root lattice -D4 described above, see |14| for more. In 
particular, one can again use the recursion ()3.24|) to calculate the coefficients /j(m) directly 
from those of the icosahedral case. The last equation also permits the determination of the 
asymptotic behaviour, compare Appendix]. With the methods described in [23 PP- 29-31], 
applied to L = Q(r), one can calculate the values of ( L (s) explicitly for s S {2,4,6}. This 
finally gives that the number of CSMs of index < is asymptotically 

|^V51og(r)C L (3)A 3 ~ 0.19773 A 3 , 

2b0 7T lz 

where Cl(3) — 1.02755 has to be calculated numerically. 

6. Concluding remarks 

In this contribution, we have shown how the so-called coincidence problem can be refor- 
mulated in a mathematical setting and then solved algebraically in dimensions 2, 3 and 4. 
Various examples have been treated explicitly, and it remains a simple exercise to work out 
tables of all coincidence rotations with small indices that could be relevant experimentally. 

Rather obvious is the question for generalizations to higher dimensions. One might hope 
that at least the root lattices could be treated in full generality, but there are complications 
from various sources. First of all, we do not have suitable generalizations of quaternions at 
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our disposal (and they proved extremely handy in our treatment), and second, we depended 
on unique factorization in one way or another - and this does not generalize to arbitrary 
lattices or modules. 

Another obvious question emerges from the observation that we have so far only dealt with 
linear isometries, while for various reasons affine extensions are necessary, in particular for a 
satisfactory formulation of the problem in the context of more general Delone sets. Though 
some preliminary investigations exists [10], more has to be done in this direction. 



Appendix A. The number of sublattices of a given index 

Given a free Abelian group of rank n, one might like to know how many different subgroups 
of (finite) index m exist. Of course, they are free and of rank n again, but here we want to 
count them separately, not up to isomorphism. Let us call that number f n (m) and derive a 
recursion relation for it 4 . Since, for fixed n, f n (m) is a multiplicative function in m, this will 
allow the derivation of a closed formula both for f n (m) and for its Dirichlet series generating 
function. 

Since any free Abelian group of rank n is isomorphic to Z n , we can treat the latter case 
without loss of generality, but with some benefit from the geometric setting. Let r be a 
sublattice of Z n of index m. Next, define a new lattice A := T D {x n = 0} by intersection 
with an (n — l)-dimensional hyperplane. Then, A is a sublattice of Z n_1 of finite index 
[Z™ -1 : A] = d, where we must have d \ m. 

At this point, we also know that r can be generated by A and some vector (y,m/d) with 
y £ Z n_1 (this is nothing but the completion theorem for bases applied to this situation). 
Next, we observe that we can actually calculate the number of lattices r that give rise to the 
same A, 

\{r | [Z n : F) = m and F n {x n = 0} = A}\ = ({choices for y (mod A)}\ = d, 

because d is the number of residue classes of A in Z n_1 . Now, summation over all possibilities 
for A results in the following simple recursion formula which is well-known, see [491 § 63, A. 13 
on p. 251], but rather difficult to locate: 

(A.l) Urn) = V d ■ f n ^(d) = m ■ V ^- • / n _i(d) . 

^-^ ^-^ m/d 

d\m d\m 

One can now derive a closed expression for f n (m), namely 
Proposition A.l. f n (m) = T,d v ..,d n =m d? • 4 ' ■ • • ' < _1 ■ 

Here, the sum runs over all n-tuples (d±, . . . , d n ) of positive integers subject to the restric- 
tion that d\ ■ . . . ■ d n = m. 



4 This derivation is partially based on notes by P. A. B. Pleasants 
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Proof: It is clear that f\{m) = 1. From Eq. (|A.1|I . we get by induction 

f n+1 (m) = Y J d -fn(d) = Y,{ d E 4-d\-...-d n n - 1 ) 

d\m d\m d\-...-d n =d 

= E E d\-dl....<r n 

d\m d\-...-d n =d 

E 4-d\-dl-...-di 

d$-d\- ...-d n =m 

which completes the argument, see [281 16U| for alternative approaches. □ 

Let us determine a generating function for f n {m)- Due to the multiplicativity of /n(m) in 
m, one would like to have a Dirichlet series generating function. This can be found as follows, 
using again the recursion relation (|A.1|) . 

„ f , n _ f /»W _ ^ Ed\m d - fn-i{d) 
n[S+ ' ~ ^ m'+ 1 ~ ^ m'+ 1 

m=l m=l 

E Y^d\m4i' fn-l{d) ^ ^ l/m ^ fn-l(£) 

m s ~ 2^1 m s 2-> £s 

m=l m=l £=1 

oo ^ 

= E^TT- F -l( S ) = C( S + l)^n-l( S ). 

771=1 

The middle line is the product formula for two Dirichlet series generating functions, applied 
to our special case. From this calculation, one gets the recursion 

(A.2) F n (s) = ((s) ■ F n ^(s - 1) . 

Since there is only one sublattice of index m for the case n = 1, we have F\ (s) = J2m=i 1 / mS = 
£(s) and thus, by induction, one obtains 

Proposition A.2. F n (s) = CO) • CO - 1) • . . . • CO - n + 1) . □ 

In particular, this gives i^O) = CORO ~~ 1); which is the well-known generating function 
for the divisor function / 2 (m) = a 1 (m) = Yld\m,d- 

Let us close this appendix with a short remark on the asymptotic behaviour of the co- 
efficients. F n (s) has its rightmost pole at s = n, with residue r 1 = 1 (if n = 1) and 
r n = C(2) • C(3) • . . . • C( n ) (if n > 1). Then, the number of sublattices with index < N 
is asymptotically given by r n ■ N n /n, while the average number of sublattices with index m 
grows like m n_1 . 
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